The Union Jack model: a quantum spin model with frustration on the square lattice 
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A new quantum spin model with frustration, the 'Union Jack' model on the square lattice, is 
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I. INTRODUCTION 

Frustrated lattice spin models in two dimensions have 
attracted much discussion in recent years. They ex- 
hibit new and interesting phase structures and phase 
transitions; in particular, they may develop 'spin liq- 
uid' states, without long-range order It is also be- 
lieved that frustrating interactions may play a role in 
the high-temperature superconducting cuprate materials. 
Primary examples are the anisotropic triangular lattice 
Heisenberg antiferromagnet (Fig. ^), the square lattice 
Ji — J2 model (Fig. QJ)), and the Shastry-Sutherland 
model (Fig. ^). In this paper we discuss another mem- 
ber of this group, the 'Union Jack' model (Fig. ^i); 
which is another square lattice Heisenberg antiferromag- 
net with frustration, and might be expected to display 
some interesting properties. 

The spin- 1/2 Ji — J2 model on the square lattice Q 
involves antiferromagnetic Heisenberg spin interactions 
with coupling Ji between nearest neighbours, and cou- 
pling J2 between diagonal next-nearest neighbours, as 
illustrated in Figure^). The Union Jack model has J2 
interactions on only half the diagonal bonds, in the pat- 
tern shown in Figure^), so that the lattice consists of 
two different site types, A and B, and the unit cell is 2x2 
sites. Both models will exhibit quantum phase transi- 
tions as the coupling ratio a — J2/J1 is varied. 

In the Ji — J2 model at small a, the Ji coupling is 
dominant and produces antiferromagnetic Neel ordering 
of the spins. At large a, the J2 interaction is dominant, 
and produces antiferromagnetic ordering on the two diag- 
onal sublattices; and then the effect of the Ji interaction 
is to align the two sublattices to form a columnar ordered 
state as illustrated in Figure [3 an example of the 'or- 
der by disorder' phenomenon. Numerical investigations 
miEllllllllSIIl have shown that the boundaries 
of these two phases lie at a ~ 0.38 and a ~ 0.60, respec- 
tively. 

The nature of the intermediate phase or phases re- 
mains controversial. Monte Carlo simulations are ham- 
pered by the 'minus sign' problem; exact diagonaliza- 
tions are limited to small lattices; and series expansions 
are based on some particular ordered reference state, and 
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FIG. 1: Lattice spin models with frustration. Solid lines rep- 
resent nearest-neighbour antiferromagnetic interactions J\\ 
dashed lines represent next-nearest-neighbour interactions J2 . 
Case: a) anisotropic triangular lattice; b) Ji — J2 model; c) 
Shastry-Sutherland model; d) Union Jack model. 



are only valid within a single phase. It is generally be- 
lieved that the intermediate phase is gapped and shows 
no long-range magnetic order. Field theory approaches 
1^ and dimer series expansions H, S H seem to 
indicate a columnar dimerized state in the intermediate 
region, with spontaneous breaking of translational sym- 
metry, as illustrated in Figure |21 Capriotti, Sorella and 
coworkers, on the other hand, have suggested a homo- 
geneous spin-liquid 'plaquette RVB' state 0, ^^id 
have found that exact diagonalization upjto 6x6 sites 
shows no strong evidence of dimerization 12]. Another 
Monte Carlo study has suggested a columnar dimer state 
with plaquette- type modulation ,13] . Sushkov et al. 
have even suggested that there may be three different 
phases in the intermediate region: reading from left to 
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FIG. 2: Phase diagram of the Ji — J2 modeL 



right, a Neel state with columnar dimerization; a colum- 
nar dimerized spin liquid, and a columnar dimerized spin 
liquid with plaquette-type modulation. 

Several discussions have centered on the Lieb-Schulz- 
Mattis theorem in higher dimensions [T^ ITEL ITfil | , which 
shows that for a spin system with half-integer spin per 
unit cell, there is an excitation energy gap behaving like 
where L is the linear size of the system. Takano et 
al. [131 have argued that a uniform RVB state without 
gapless singlet excitations is excluded by the theorem, 
and that the true ground state is a plaquette state with 
spontaneously broken translation invariance and fourfold 
degeneracy. Later arguments 0, 0, 0| have refuted 
this, however, and shown that the theorem may be sat- 
isfied if the translation symmetry remains unbroken, but 
the ground state has a fourfold topological degeneracy 
instead, as in a simple dimer model. 

The Hamiltonian for the Union Jack model is 



H 



(1) 



A: <7i7in> 



where the Ji and J2 interactions connect sites as shown 
in Figure ^i)- A classical variational analysis shows that 
for a = J2/J1 < 0.5, the ground state is the Neel state, 
as in the Ji — J2 model. For a > 0.5, the ground state is 
the canted ferrimagnetic state shown in Figure 13 where 
the spins on the A sublattice are canted at an angle 9 to 
those on the B sublattice, and 26 to their neighbours on 
the A sublattice. The energy of this state is 



Eo = NS^ [a cos 261-2 cos 9] 



(2) 



where we have set Ji = 1, 5 is the total spin per site, 
and N is the number of sites. This energy is minimized 
when 



sin6'(2acose' - 1) = 0. 



(3) 



For a < 0.5, the lowest energy corresponds to sin 6* — 0, 
i.e. the simple Neel state. For a > 0.5, the lowest energy 
solution is 



cos 9 — —, 
2a 



(4) 



corresponding to the canted state. In the limit a 00, 
the angle 9 — > 7r/2: the spins on the A sublattice are 
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FIG. 3: Spin configurations in the canted phase. 



Neel ordered, as expected, and the spins on the A and B 
sublattices are at right angles to each other. 

In the canted phase, there is a net magnetization per 
site in the z direction, 



S 



S . 



M, = -(l-cos0) = -(l- 
and a staggered magnetization per site 



2a' 



S 

Ms = -(1 + cos( 



^(1 + ^). 
2^ 2a' 



(5) 



(6) 



The model thus demonstrates a new mechanism for 
causing ferrimagnetism. Ferrimagnetism usually arises 
when the individual ionic spins are different on differ- 
ent sublattices. Here the spins have the same magni- 
tude on each sublattice, and the interactions are antifer- 
romagnetic, but the frustration between them produces 
a nonzero overall magnetic moment. 

In the neighbourhood of the transition point, ac — 0.5, 



we find for the ground state energy per bond eo 
in zero magnetic field : 



3^ 

' 8 ' 



da 



d^eo 



0, a ^ ac— 

y2 



Eo/2N 

(7) 
(8) 
(9) 



da^ { -45-^, a ac+ 

while for the net magnetization in the z direction, 

M^r^ S{a-ac), a->ac. (10) 

Thus the transition is not quite the simple first-order 
transition found classically for other frustrated models. 
There is no "latent heat", i.e. no discontinuity in deo/da, 
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although there is a finite discontinuity in the second 
derivative. The magnetization vanishes hnearly at 
ac, and the energy gap is zero on both sides of the transi- 
tion, corresponding to spontaneous breaking of the spin 
rotation symmetry Translation symmetry on the lat- 
tice is also broken on both sides of the transition. The 
question then arises, whether quantum fluctuations will 
change the character of the transition. 

In the remainder of this paper, we shall discuss how 
this picture is modified when one goes beyond the clas- 
sical analysis to a spin- wave treatment. In Section UTI we 
shall present a modified spin-wave treatment of the Neel 
phase to second order; and in Section UTTl we shall give 
a spin-wave treatment of the canted phase to first order. 
The results are discussed in Section Hvl Our conclusions 
are summarized in Section Ivl 



|20l l2ll | : this requirement turns out to be equivalent to 
Takahashi's modified spin- wave theory . A very sim- 
ilar treatment can be applied to the Union Jack model. 
We will follow Gochev's notation as far as possible. 

Let us rewrite the Hamiltonian f^l, adding in a stag- 
gered magnetic field h: 

<lm> <lil2> I ™ 

(12) 

where we have divided the lattice into even (A) and odd 
(B) sublattices, denoted by indices I and m respectively, 
and set Ji = 1. 

Introduce creation and annihilation operators for the 
"spin deviation" on the two sublattices by means of a 
Dyson-Maleev transformation: 



II. SPIN- WAVE THEORY FOR THE NEEL 
PHASE. 

A. Formulation 

A review of spin-wave approaches to the Ji — J2 model 
has been given by Gochev 19] . Chandra and Doucot 
gave the first conventional spin- wave treatment. The 
results of this approach were found to become unstable 
as the frustration a increases, due to strong interactions 
between the spin wave bosons. The interactions involve 
a quadratic term 



Ho = 



k 



(11) 



which is comparable with the zeroth-order Hamiltonian 
for 5 ~ 1 and a ~ 1. Thus one needs to go to a mod- 
ified treatment in which this quadratic term is absent 



Sf 



S — a]ai 



S+ = {2Sy/^ai-{2S)-'/^alaiai 
and then perform a Fourier transform 



(13) 



Ok 



1/2 



ikl 

e ai; 



1/2 



-ikm I 



(14) 



to give the Hamiltonian in the form 



H 



2S 



)+4S{{l-a 



k k k 



— Y '^i-2-3+4[273-4ala24^4 + 74ala2a3fe4 + 7iai44^4] 



^ Y ^i+2-3-444a3a4('72~4 - ^{m + V2)) 



(15) 



r 



where we have used the shorthand notation 1 • • • 4 for 
momenta ki • • •k4, and 7k,?7k are the structure factors 
for the full lattice and the A sublattice respectively: 



7k 



COS — cos -r 
2 2 



-(cosfca; -I- COS ky), 



(16) 



if we set the spacing of each sublattice equal to 1, and 
kx , ky are the components of momentum along the diag- 
onal axes of the two sublattices. 

The Hamiltonian can now be diagonalized up to second 
order by a Bogoliubov transformation: 

Ok = Wkak - -Wk/Sk, ^k = "k/?k - V]^al. (17) 
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where 



2 2 1 



(18) 



After normal ordering the transformed Hamiltonian, 
the condition that off-diagonal quadratic terms vanish 
turns out to be: 

Qk = 45[7k(wk + "k)-2TO(l-^-|(l-?7k))] 
+4(i?i - R2)[jk{ul + vl) - 2ukt'k] 
-4auki;k(i?2 - -R3)(l - Vk) 







(19) 



where the lattice sums Ri are 

2 

N 



2 2 

-Rl = T7 XI ^kUkWk, ^2 = ^ 51 ^3 



TV 



(20) 

A solution to equation (|19|l can easily be found (for 
h < 0): 



Uk = 



l+Ek 



2ek 



1/2 



, Vk = sgn(7k) 



1-ek 



-,1/2 



2ek 



(21) 



P = 



5 + i?3 — i?2 
5 + i?i — i?2 



(23) 



(the parameter p is Gochev's a). Then one merely has 
to find self-consistent solutions for the two parameters a 
and p, given by equations l|23|) . 

We can now write the Hamiltonian as 



Hdm — Wq + Ho + Vdm 
The constant term is 
Wo 



(24) 



27Veo = 2N[-{S + Ri - R2) 

fy J-f 

-{S + Rs - R^f + -{S - R2)] (25) 



(note that there is a misprint in Gochev [T3 at this 
point). The quadratic part Hq is diagonal: 



^o = X(^k4«k + i?f/3i/3k) 



(26) 



where the a and /3 bosons have different spin- wave ener- 
gies in this case: 



where 



Ek 
/k 



= 1- 



I -ha- p-{l - 77k) 



1 

i{S + Ri-R2) 



(22) 



= ^S[ulil-ail-r,k)-—) + vlil-—) 

-2ukWk7k] + 4(i?i - i?2)(Mk + "^k ^ 2-UkWk7k) 

+Aaul{R2 - RsKl - Vk) (27) 
and is a similar expression with and ?;k inter- 
changed. 

The normal-ordered quartic interaction operator Vdm 



Vdm - X ^1+2-3-4 [$^'^aia2/33/?4 + ^^^^ aialpllSl - 2$(3)at ^^^.^^.^ _ 2$(4)^t ^t^t^g 

1234 

-2$(5)/3t^3^^Q,^ _ 2$(6)a^44a2 + «'^^^/?I4/33/34 + ^'^^^ alala^a2 + 4$(9)4/3| 02/33]- (28) 

I 



Explicit expressions for the vertex functions $(*)(1234) 
of the Neel phase are given in Appendix A. 



B. Higher-order corrections 



We can now use perturbation theory to calculate the 
next-order (1/5*) corrections to these results. The lead- 
ing correction to the ground-state energy corresponds to 
Figure^), and is given by 



To order 1 in a l/S" expansion, the ground-state energy 
per bond is eg, as given by equation The staggered 
magnetization per site is 



M, = 2^\h=n = S - R2, (29) 
and the spin- wave energy is that given by equation (|27|l . 



AEq 
2N 



1+2-3-4- 



4>(i)(1234)4>(2)(34i2) 



1234 



(30) 

while the corrections to the spin-wave energies E^ and 

4 



Ek (Fig. W) are 
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Ai?" = _8(-VV'52 3 1 ^ 'i^^'Hlk23)$(3)(23kl) J2\'^^ ^ ^_<i>W(123k)$W(k321) 

^ \ N / ^-^ ■^^^ I- P/3 I T^Q I 771(1 77Q\ \ N ] ^-^ ^^■^ 



R/'^y^ $(^)(lk23)<i>W(23kl) /2y^ $(i)(i23k)<i>(^)(k321) 




(a) (b) 



FIG. 4: Perturbation diagrams in the Neel phase. Solid lines 
correspond to a bosons, dashed lines to /3 bosons, a) Leading 
correction to the ground-state energy, b) Leading corrections 
to the spin- wave energy 



The corrections to the staggered magnetization can 
also be represented in terms of perturbation diagrams 
as follows. We have 



= ^ < V^ol E - E > (32) 

; in 

= < Vo|5-i?2-V^M-V//|V'o > (33) 

where 



yli = E("k + ^'k)K"k + /3^/?k) (34) 

k 

= E^^k^'kCai./Sl + ak/Sk) (35) 

k 

Hence the leading corrections to of 0(1/5") corre- 
spond to the diagrams shown in Figure 13 and are given 

by 



AM,, 



16 ^ $(i)(1234)$(2)(34i2) 



1234 



1111 
ei 62 £3 £4 



F(1234) 



32_ ^ 



(36) 
(37) 



where 



F(1234) = <I>(i)(1234) 



7i$('^)(341 2) 73$(6)(4312) 



fiei{E^ + E^) he^{E% + Ei 



$(2^(3412) 



7i$(5)(2i34) 73$(3)(1234) 



he^{E1+E^^) faeaiES + E^ 



r 



m 



The results of these calculations are discussed in Sec- to be somewhat more complicated in this case, and we 
tion llVI shall only keep terms down to 0(S) in the energy. We 

write the Hamiltonian as 



III. SPIN- WAVE THEORY FOR THE CANTED 
PHASE. 

In this section we shall present a modified spin wave 
treatment of the canted phase. The formalism turns out 



H 



E Si-Sj+a E Si-Sj 



A: <nnn> 



-/iiE^f + /^2E^^ 



(39) 
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quantizing the spins with respect to the configuration 
shown in Figure 3, so that the z axis on the B sublattice 
points upwards, and the z axes on the sublattices Ai and 
A2 are canted at angle Q to the downwards direction as 
shown, where is a parameter to be determined. The y 
axes are taken to he perpendicular to the paper in each 
case. In terms of spin components, we then have: 




FIG. 5: Perturbation diagrams contributing to the Neel mag- 
netization at order l/S. 



^ [('5'll:n'5'l2:n+;i' + '5'll:n'^l2:n+Ai' ) COS 20 + '5'^l:n'S'A2:n+Ai' + {'^\\:J^\l:n+ li.' ~ '^\\:J^\l:rx+ li.') "^"^^ ^0 
Al:n,/i' 

+ ^lE^B:n + /l2 5]^in (40) 
B:n A-.-a 



r 



where the direction vectors are {/x} = ±i, ±j; {/x'} = 
±(i ± j) and the phase factor r\■n^x = ±(—1)"" for /i equals 
±i or ±j, respectively. 

Now we introduce boson creation and annihilation op- 
erators for the spin deviation on the two sublattices by 
means of the Dyson-Maleev transformation: 





— 5* ojj^ajn 






'-'Aim 












-{2S) 


"m"m"m 




= s-bibn 










{2S)- 





^B:n = C^Sy/'bl (41) 

and perform a Fourier transform 

\ / 

/ o \ 1/2 

(n) E^*'"^- (42) 

to obtain 
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H 



NS^ 2 cos 6* -a cos 261 



25* 



^ + 2(7V53)i/2(sin6l - a sin 261) (a^,^ + a^^^) 



-2S[{2 cos - ^) ^ 6t &k + (2(cos - a cos 20) - ^) ^ aj, Ok - (1 + cos 0) ^ 7k(&ka-k + ?>l,a^_k) 

k k k 

-(1 - COS e) ^ 7k(&kak + ^k^k) + ^ XI '^k[2(l + cos26l)aJ^ak - (1 - cos 26l)(aka_k + ala^_^)]] (43) 
k k 

I 



to order 5, where as usual 



7k 



'7k = 



,ik-p' 



(44) 



The term of order S"^ is simply the classical energy of 
this configuration, equation The term of order S*^/^, 
proportional to (fo^^Tr + ^t ■^)^ would indicate that we have 
not chosen the optimum reference state, and should be 
set to zero: this gives a condition on the angle identical 
to the classical condition, equation At higher orders, 
this condition would give a 'renormalised' value for the 
angle 9. 

The term of order S is quadratic in fermion operators, 
and can be diagonalized by a Bogoliubov transformation. 
In the sector involving momenta (±k), the quadratic part 
of the Hamiltonian is 



i/k = 2S[{ al fe_k hi a_k ) ^k 



/ Ok \ 
&k 

V «^-k J 



+ iVk] (45) 



where to leading order the matrix iJk has elements hij 
given by 

^2 

hii = /144 = 2(cos — a cos 26*) + a(l + cos 26')77k — — 



^22 = ^33 = 2*^°^^ ~ 



25' 
(46) 



(47) 



hi2 = = /134 = ^43 = -(1 + cos6')7k; (48) 
hi3 = hzi = /i24 = /i42 = (1 - cos6')7k; (49) 
h2z = h^2 = 0. (50) 

hii = /i4i = — a(l — cos20)77k; (51) 

Note that the matrix is symmetric about both diago- 
nals. The normal-ordering correction to 0{S) is 



iVk 



-/ill 



(52) 



Thus the procedure here involves a 4x4 matrix diago- 
nalization, rather than 2x2. We want a transformation 
which preserves commutation relations, and the symme- 
tries of the problem, i.e. we need to find a symplectic 
transformation: 



(53) 



/ ak \ 




/ "k \ 


&ik 
6k 


= 5k 


/?lk 
/3k 


V «'-k ) 




Uik/ 



where the elements of Sk obey 



z^4 



3 



^ 



Then 



S HkS. 



(54) 
(55) 



There are initially 16 unknown parameters, corre- 
sponding to the elements of the transformation matrix 
^k. The conditions H54() turn out to eliminate 12 of these, 
leaving only 4 independent parameters; these 4 parame- 
ters are determined by the condition that H^^ should be 
diagonal, i.e. 



^12 ^ ^13 ^ ^14 ^ ^23 ^ 



(56) 



We have determined the solution to this problem numer- 
ically. In its diagonalized form, the Hamiltonian now 
reads: 



+h'22WlPi. + PUP-k) + K], (57) 
Ni, = /I'li + - /ill - h22 (58) 



where 



Hence one can easily obtain numerical results for the 
ground state energy and single-particle dispersion rela- 
tions. The staggered magnetization is given by 



Ms - ^<^0\J2^B 
B 



cos ( 



A 



TABLE I: Ground-state energy per bond as a function of 
a. Linear spin-wave theory: Eq^'; modified second-order 
spin- wave theory: Eq^' ; with higher-order corrections eg''' . 



Q 







Neel expansion 




^(3) 



Canted expansion 


0.0 


-0 


32897 


-0.33521 


-0 


33503 




0.1 


-0 


31778 


-0.32523 


-0 


32496 




0.2 


-0 


30666 


-0.31543 


-0 


31509 




0.3 


-0 


29561 


-0.30586 


-0 


30543 




0.4 


-0 


28465 


-0.29652 


-0 


29602 




0.5 


-0 


27377 


-0.28746 


-0 


28691 


-0.27377 


0.6 


-0 


26300 


-0.27870 


-0 


27817 


-0.25465 


0.7 


-0 


25233 


-0.27028 


-0 


26990 


-0.24458 


0.8 


-0 


24180 


-0.26224 


-0 


26224 


-0.24050 


0.9 


-0 


23141 


-0.25462 






-0.24058 


1.0 


-0 


22119 


-0.24746 






-0.24367 
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FIG. 6: Ground-state energy per bond as a function of cou- 
phng Q. Short-dashed hue: classical result; long-dashed line: 
linear spin-wave theory; solid line: 2nd order spin-wave the- 
ory. 



TABLE IL Net magnetization M-^ and staggered magnetiza- 
tion Ms as functions of a. Conventions as for Table |T] 



a 






Neel expansion 

Mf^ 






Canted 


expansion 

mV 











30340 


0.30340 





30601 









1 





29582 


0.29149 





29499 









2 





28759 


0.27819 





28265 









3 





27859 


0.26327 





26872 









4 





26873 


0.24647 





25270 









5 





25785 


0.22751 





23386 


0.25785 


0.00000 





6 





24578 


0.20607 





21088 


0.28800 


0.02688 





7 





23230 


0.18181 





18163 


0.30067 


0.05143 





8 





21714 


0.15442 





14529 


0.30525 


0.07207 





9 





19992 


0.12369 






0.30579 


0.08912 


1 








18017 


0.08954 






0.30429 


0.10321 



5* 



(1 + cos( 



k 



1-f 4+cos0(5?2 + s?4)] (59) 



and the net magnetization in the z direction: 



N 

|(1-COS0) 

2 |. 2 2 

~Jf + ^23 

k 



cose{Sf^ + sli)] (60) 



IV. NUMERICAL RESULTS 

Numerical results for the model have been obtained 
using the finite-lattice method. The momentum sums 
were carried out for a fixed sublattice size L — y/N/2, 
using discrete values for the momenta and ky, e.g. 



kx{i) 



27r(i - 1/2) 



i = 1, • • -i, 



(61) 



where we take half-integer values corresponding to anti- 
periodic boundary conditions to avoid any integrable 
singularities at k = 0. Results were obtained for L = 
8, 9 • • • 12, and a fit in powers of 1/L was made to extrap- 
olate to the bulk limit L ^ oo. The finite-size corrections 
for the ground-state energy per bond scale asymptoti- 
call y li ke and those for the magnetization like 1/L 

|24Ll25l |. The resulting bulk estimates are shown in Tables 
U and |n] The values at a = agree with those obtained 
by Gochcv 19J for the pure Heisenberg case. 

Figure shows the behaviour of the ground-state en- 
ergy resulting from the classical calculation and linear 
spin-wave theory for both the Neel phase and the canted 
phase, and the modified treatment with corrections for 
the Neel phase. In linear spin-wave theory, the Neel and 
canted results coincide at a = 0.5, as they must do since 
9 = there. Somewhat surprisingly, however, the Neel 
energy remains the lower of the two beyond that point, 
until the two curves cross once more at a ~ 0.84. In 
other words, the transition between the Neel and canted 
phases is pushed out to a ~ 0.84 in linear spin- wave the- 
ory, and is clearly first-order. The modified treatment 
to second order lowers the energy a little further, while 
the higher-order corrections are virtually negligible, and 
indistinguishable on the diagram. 

FigureHshows the magnetizations Ms and Mz as func- 
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tions of the coupling a. In linear spin-wave theory, the 
staggered magnetization Alg in the Neel phase is reduced 
by quantum fluctuations, as expected, but the effective 
coupling is only half that in the Ji — J2 model, and so Mg 
remains substantial at a = 0.5. In the modified second- 
order treatment it is lowered somewhat further, while the 
higher-order corrections are small and positive for low a, 
and turn negative beyond a 2± 0.7. The staggered magne- 
tization in the canted phase is also reduced by quantum 
fluctuations, although the effect is reduced at large a. 

Figure |H1 illustrates the spin- wave dispersion of the a 
and f3 bosons in the Neel phase as given by the second- 
order theory. It can be seen that the dispersion curve for 
the p bosons remains virtually unchanged at all a. That 
for the a bosons, however, develops an instability at k = 
(tt, tt), and the energy gap is predicted to vanish when 
a gets too large, signalling a transition. The instability 
occurs at a = 1.0 in linear spin wave theory, and a = 
0.645 in the modified second-order theory. 

Figure El shows similar plots from the linear theory in 
the canted phase. The energy gap for the /? bosons re- 
mains zero at all couplings, with linear dispersion near 
the origin; but as soon as we move away from a = 0.5 
(i.e ^ = 0) an energy gap opens up for the a bosons, 
with quadratic dispersion near the origin. This is in ac- 
cordance with the general counting of Goldstone modes 
in non-relativistic systems ,26] , and with previous discus- 
sions of Goldstone modes in canted phases ^3, HE @| , 
where the spontaneous symmetry breaking pattern is 
SU(2) Z2. 



V. CONCLUSIONS 

The Union Jack model is another example of a spin- 
1/2 Hcisenbcrg antifcrromagnet with frustration on the 
square lattice. A classical variational analysis predicts a 
transition at a = 0.5 between the Neel phase at small 
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FIG. 8: Dispersion curves in the Neel phase along high- 
symmetry lines in the Brillouin zone for various a. Curves 
a) and b) are from linear spin-wave theory, c) and d) from 
the modified second-order theory. Cases a) and c): a bosons; 
cases b) and d): /3 bosons. 
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a, and a canted ferrimagnetic phase at large a. Thus 
the system exhibits a new phenomenon, namely fcrri- 
magnetism induced by frustration. 

In linear spin-wave theory, the transition is pushed out 
to Qfc — 0.84. This makes the model potentially much 
more interesting, because it is very likely that the Neel 
magnetization will vanish at or before that point. In the 
Ji - J2 model, numerical studies [111. Isl. ItL H H ITolIT^ 
show that the Neel magnetization vanishes at about 
Oc — 0.38. The frustration in the Union Jack model 
is only half as strong, as seen by a glance at Figure 
or from linear spin-wave theory, and so one might ex- 
pect the magnetization to vanish at about ac — 0.76 in 
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this model. A plausible scenario, then, is that a second- 
order Nccl transition might occur at ~ 0.76, possibly 
followed by an intermediate spin-liquid phase as in the 
Ji — J2 model, and then a first-order transition to the 
canted phase at somewhat larger a. Numerical experi- 
ments would be necessary to ascertain if this is indeed 
the case. 
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APPENDIX 

The vertex functions i — 1 • • ■ 9 are: 

(1234) = 7(4 — l)viU2V3U4 + 7(4 — 2)uiV2V3U4 + 7(3 — l)viU2U3Vi + 7(8 — 2)uiV2U3V4 

-7(4)uiU2l'3W4 - ^{3)uiU2U3Vi - ^{3)viV2V3Ui — ^{4:)viV2U3V4, (62) 

—aQuiU2V3V4, 

$(2) (1234) = 7(4 — 2)viU2U3V4 + 7(3 — 2)viU2V3U4 + 7(4 — l)uiV2U3V4 + 7(8 — l)UiV2V3U4 

-7(3)^1 f2W3W4 - 7(4)^1 - 7(3)miU2M3V4 — ^ {4:)uiU2V3U4 (63) 

—aQviV2U3U4 

$(3) (1234) = 7(4 — l)viU2U3U4 + 7(3 — l)ui'U2^^3i'4 + 7(4 — 2)uiV2U3U4 + 7(3 — 2)'Uit;2t;3t;4 

-7(4)uiU2'"3'W4 - 7(3)uiU2t^3^^4 - ^{S)viV2U3U4 — 'y{A)v\V2V3V4 (64) 

— aQuiU2U3V4 

$("')(1234) = 7(3 — l)uiV2U3U4 + 7(4 — 1)miZ)2W3W4 -|- 7(3 — 2)t>iU2U3M4 -|- 7(4 — 2,)viU2V3V4 

-^{2,)ViV2U3U4 — 'y{A)v\V2V3V4 - ^{2,)uiU2V3V4 - 7(4)uiW2«3W4 (65) 
—aQv\V2V3U4 

$(5) (1234) = 7(4 — l)viV2V3U4 -I- 7(4 — 2)uiU2V3U4 + 7(3 — l)viV2U3V4 -\- 7(3 — 2)ui'U2W3t;4 

-7(4)uiU2i'3ii4 - 7(3)uiti2U3^^4 - 7(3)t;iU2t;3W4 - 7(4)tiiW2W3i^4 (66) 

-aQuiV2V3V4 

$(6) (1234) = 7(4 — 1)U4U2U3V4 + 7(4 — 2)ui'U2W3t;4 -|- 7(3 — l)uiU2V3U4 + 7(3 — 2)?;it;2t;3U4 

—^{^)viU2U3V4 — ^{2,)V\U2V3U4 — ^{2,)u\V2U3V4 — 'y{A)uiV2V3U4 (67) 

— aQviU2U3U4 

$^"^^(1234) = 7(4 — 2)viU2V3U4 + 7(4 — l)uiV2V3U4 + 7(3 — 2)viU2U3V4 + 7(3 — l)uiV2U3V4 

■~^(4)V\V2V3U4 — j{3)ViV2U3V4 - j{3)uiU2V3U4 — 7(4)uiM2W3t'4 (68) 
-aQviV2V3V4 

$(^) (1234) = 7(4 — 2)uiV2U3V4 + 7(4 — l)viU2U3V4 + 7(3 — 2)uiV2V3U4 + 7(3 — l)viU2V3U4 

-j{'i)uiU2U3V4 — ^{3)uiU2V3U4 - 7(3)t;iW2W3l^4 " ^{4)viV2V3U4 (69) 

—aQuiU2U3U4 

$(9) (1234) = 7(4 - 2)uiU2U3U4 + 7(4 — l)viV2U3U4 + 7(8 - 2)uiU2V3V4 + 7(8 - l)wit'2f3t'4 

— 7(4)uiU2W3U4 - 7(3)miI'2U3^'4 — ^{3)viU2U3U4 - 7(4)t)iU2V3l'4 (70) 
-aQviU2V3U4 

where 



Q = 7?(3-2) + r,(4-2)-r,(3)-r,(4). 



(71) 
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